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The High Altitude Water Cherenkov (HAWC) gamma-ray observatory 
is a wide £eld-of-view observatory sensitive to 100 GeV — 100 TeV gamma 
rays and cosmic rays. The HAWC observatory is also sensitive to diverse 
indirect searches for dark matter annihilation, including annihilation from 
extended dark matter sources, the diffuse gamma-ray emission from dark 
matter annihilation, and gamma-ray emission from non-luminous dark 
matter subhalos. Among the most promising classes of objects for the 
indirect detection of dark matter are dwarf spheroidal galaxies. These 
objects are expected to have few astrophysical sources of gamma rays but 
high dark matter content, making them ideal candidates for an indirect 
dark matter detection with gamma rays. Here we present independent 
limits on the annihilation cross section for 14 dwarf spheroidal galaxies 
within the HAWC £eld-of-view, as well as their combined limit. These 
are the hrst limits on the annihilation cross section using data collected 
with HAWC. 
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1 Introduction 


Evidence for the existence of dark matter has been seen in many observations, such 
as in galactic rotation curves, galaxy clusters, gravitational lensing, large-scale cos¬ 
mological structure, and the cosmic microwave background. While there is ample 
evidence of the gravitational effects of dark matter, the particle nature of dark mat¬ 
ter remains unclear. Weakly Interacting Massive Particles (WIMPs) are among the 
leading hypothetical candidates for dark matter. WIMPs can annihilate into stan¬ 
dard model particles and produce photons via pion decay, radiative processes by 
charged leptons, or direct production of gamma rays through loop-order processes 
|1]. Indirect detection of dark matter can be done by looking for the products of 
dark matter annihilation in objects with high dark matter content. When there is 
no signihcant detection of these by-products we can set limits on the annihilation 
cross section. Dwarf spheroidal galaxies are prime targets for indirect searches of 
dark matter as these objects are expected to have very few astrophysical sources of 
gamma rays, but high dark matter content. Benefiting from the low signal-to-noise 
ratio for these sources, we can constrain the dark matter annihilation cross section 
for several annihilation channels. 


2 The HAWC Observatory 

The High Altitude Water Cherenkov (HAWC) gamma-ray observatory is a wide field- 
of-view observatory located at an elevation of 4100 m on the volcano Sierra Negra 
in the state of Puebla, Mexico. It is a second-generation water Cherenkov detector, 
building on the success of the Milagro experiment, and is sensitive to 100 GeV- 
— 100 TeV gamma rays and cosmic rays. HAWC consists of a densely packed array of 
300 Water Cherenkov Detectors (WCDs), each of which is 7.3 m in diameter and 4.5 m 
deep, with the total array covering an area of approximately 20,000 m^. Each WCD 
contains a light-tight hermetically sealed plastic bag, called a bladder, containing 
200,000 L of ultra-purified water and is instrumented with four photomultiplier tubes 
(PMTs). Three of the PMTs are 8-inch Hamamatsu R5912 PMTs, re-used from 
the Milagro experiment, with the fourth and central PMT being a 10-inch R7081- 
MOD high-quantum efficient PMT. The PMTs detect Cherenkov light from energetic 
particles passing through the WCD, which are produced in the extensive air shower 
of the primary incident particle. HAWC has an instantaneous £eld-of-view of 2 sr, 
and a duty cycle >95%. 
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3 Gamma-Ray Emission from Annihilating Dark 
Matter 

3.1 Dark matter differential flux 

A calculation of the expected gamma-ray flux from dark matter annihilation requires 
information about both the astrophysical properties of the potential dark matter 
source and the particle properties of the initial and hnal-state particles. The differ¬ 
ential gamma-ray flux integrated over the solid angle of the source is: 

dF _ {aAv) dN^ 

dE ~ StiMI dE ^ 

where {aAv) is the velocity-weighted cross section, dN^/dE is the gamma-ray spec¬ 
trum per dark matter annihilation, and is the dark matter mass [1]. We dehne 
the dark matter J-factor (J) as the mass density (p) squared integrated along the 
line-of-sight x, averaged over the solid angle of the observation region: 

J= f dfl f dxp‘^{rgai{6,x)) (2) 

J source J 

where the distance from the source is given by: 

rgai{0,x) = \Jr^ - 2xRcos{9) + x"^ (3) 

where R is the distance to the center of the source, and 9 is the angle between the 
source and the line-of-sight x [1]. 

3.2 Calculation of dark matter spectra 

Pythia 6.4 [13] was used in this analysis to calculate the expected photon spectrum for 
a particular WIMP annihilation channel. The photon radiation of charged particles 
was simulated, as well as the decay of particles such as the 7r° [U [Tl|. For each 
annihilation channel and each dark matter mass, the average number of photons in 
each energy bin per annihilation event was calculated. This differential flux, dN^/dE, 
was used to determine the dark matter J-factor of the targeted source. The dark 
matter annihilation channels considered in this analysis are: XX bb, XX , 

XX -t XX -t tt and yy W+W-. 

3.3 Dark matter distributious 

In this analysis there are two dark matter density prohles that are used, the NFW 
(Navarro-Frenk-White) model and the Einasto model. The density prohles describe 
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how the density (p) of a spherical system varies with distance (r) from its center. The 
NFW prohle as developed by Julio Navarro, Carlos Frenk and Simon White , is 
the simplest single parameter model consistent with N-body simulations. The NFW 
density prohle takes the form: 


PNFw{f') 


Ps 

(r/r^)(l + t/tsY 


(4) 


where ps is the scale density and is the scale radius of the galaxy. The Einasto 
prohle, as developed by Jaan Einasto [2] is of the form: 


PEinasto('^) 


PsGXp 




(5) 


where a controls the degree of curvature of the prohle. The Einasto model is a 
more conservative estimate of the dark matter density prohle than the NFW model. 
Depending on the chosen density prohle model, we can substitute p(r) into equation 
l^to calculate J for a particular source. We use the Einasto model with a = 0.303 for 
Seguel [T], while we use the NFW model for the remaining dwarf spheroidal galaxies. 
The source parameter values for the fourteen dwarf spheroidal galaxies presented in 
this analysis are listed in table [T] 


4 Analysis 


4.1 Likelihood analysis 

In order to analyze a particular region of the sky, we do a likelihood ratio test. This 
allows us to estimate the signihcance of a source that has a low signal-to-noise ratio. 
The likelihood ratio test is a ratio of two diherent hypotheses: 1) the null hypothesis, 
which assumes that no source exists and all observed photons are due to background, 
and 2) the alternative hypothesis, in which a source exists. According to [HI [12], 
we can directly relate the signihcance of the observed result to the likelihood ratio 
Lq/T, where Lq is the null hypothesis likelihood and L is the alternative hypothesis 
likelihood: 

( 6 ) 

where TS is the Test Statistic. Lq and L are maximized with respect to any free 
parameters. In this analysis we use data with no free parameters for Lq and a single 
free parameter for L. Since the null and alternative hypothesis likelihoods are taken 
to be Poisson distributions with parameter /r, then: 



TS 



(7) 
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Source 

RA 

Dec 

Ps 

rs 

R 


J 

a 

Ref. 




GeV /cvcA 

kpc 

kpc 

GeV^cm ®sr 



Bootes 1 

210.05 

14.49 

8.12 

0.27 

66 

3.8 

X 

10^^ 

-0.04 

PC] 

Canes Venatici I 

202.04 

33.57 

0.79 

0.55 

218 

2.9 

X 

10^6 

0.91 

PE] 

Canes Venatici II 

194.29 

34.32 

4.77 

0.13 

160 

2.5 

X 

10^6 

0.34 

PP 

Coma Berenices 

186.74 

23.90 

9.76 

0.16 

44 

2.6 

X 

lO^s 

0.88 

PP 

Draco 

260.05 

57.07 

0.98 

2.1 

76 

2.0 

X 

1019 

0.30 

mm 

Hercules 

247.72 

12.75 

0.80 

0.32 

132 

1.6 

X 

10^6 

-1.67 

pp 

Leo I 

152.11 

12.29 

16.20 

0.28 

254 

1.2 

X 

lO^s 

0.13 

Pi 

Leo II 

168.34 

22.13 

162.01 

0.06 

233 

1.2 

X 

ioi« 

-0.02 

pp 

Leo IV 

173.21 

-0.53 

1.99 

0.15 

154 

7.3 

X 

10^5 

0.51 

pp 

Segue 1 

151.75 

16.06 

4.18 

0.15 

23 

1.8 

X 

10l9 

-0.33 

pp 

Sextans 

153.28 

-1.59 

3.38 

0.37 

86 

1.0 

X 

lO^s 

-1.55 

pp 

Ursa Major I 

158.72 

51.94 

2.39 

0.31 

97 

2.3 

X 

1017 

-0.37 

pp 

Ursa Major II 

132.77 

63.11 

13.79 

0.17 

32 

1.1 

X 

1019 

0.10 

pp 

Ursa Minor 

227.24 

67.24 

3.89 

0.65 

76 

9.6 

X 

1018 

0.26 

pp 


Table 1: Astrophysical parameters for the fourteen dwarf spheroidal galaxies within 
the HAWC £eld-of-view and their references. The source, right ascension {RA), 
declination {Dec), scale density {ps) in GeV/cm^, scale radius (r*) in kpc, distance 
to the source {R) in kpc, and the dark matter J-factor in GeV^cm“®sr are all listed 
above. The signihcance (a) is also shown for each source. The signihcances listed are 
for = 10 TeV and the XX bb annihilation channel. Segue 1 uses an Einasto 
dark matter density prohle with a = 0.303, while the remaining sources use an NEW 
prohle. 
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For the alternative hypothesis, /i is the number of expected counts in each bin 
= E + B where E is the expected number of signal counts and B is the num¬ 
ber of background counts), and k is the number of total events in each bin from data 
{N). For the null hypothesis, we make similar definitions, where po is fhe number 
of expected counts in each bin for the null model (i?), and k is the number of total 
events in each bin from data {N). So equation becomes: 


ts = Y. 

bins 


2iV In ( 1 + ^ ) - 2E 


( 8 ) 


where TS is summed over all energy-proxy bins and spatial bins [5] . In this analysis 
the dwarf spheroidal galaxies were treated as point sources. The formalism of the 
likelihood ratio test has been implemented in the HAWC software, and is known as 
the Likelihood Fitting Framework (LiFF). The details of LiFF are described in [TO] . 


4.2 95% confidence level limit 

We use the likelihood ratio test and definition of TS* in equation to constrain 
the dark matter annihilation cross section for a particular source. To set a 95% 
Confidence Level (CL) limit, we first maximize TS {TSmax) and then optimize ATS = 
TSmax - TSq5 such that: 
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TS 


max 


E 

bins 


2Nln 



— 2^Eref 


(9) 


where ^ is the scale parameter with which we scale our expected counts from a source 
until we are 95% confident that HAWC will be sensitive to that particular dark matter 
annihilation cross section. Eref is the expected number of counts in each bin for the 
reference cross section, ((TA'y)ref- The limit is not dependent on the value chosen for 
{aAv)ref as our expectation is linearly proportional to the cross section, as can be 
seen in equation The scale parameter ^ is then used to set the 95% CL limit for a 
particular source, dark matter annihilation channel and dark matter mass: 


{(^Av)95% = ^ X {aAv)rei 


( 10 ) 


5 Dark Matter Limits from HAWC 

We present in this analysis individual limits from the fourteen dwarf spheroidal galax¬ 
ies within the HAWC £eld-of-view. The limits shown here are done assuming the 
dwarf spheroidal galaxies act like point sources, not extended sources. Through 
detailed simulation of the HAWC gamma-ray sensitivity and backgrounds [5], we 
have determined the signihcance of the gamma-ray flux for a range of dark matter 
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Figure 1: The dark matter annihilation cross section limits for the fourteen dwarf 
spheroidal galaxies within the HAWC £eld-of-view. Individual limits are shown for 
two select dark matter annihilation channels: XX bb and XX The solid 

black line is the combined limit from a stacked analysis of all fourteen dwarf spheroidal 
galaxies for the shown annihilation channels. 


masses, 0.5 TeV — 1000 TeV, and five annihilation channels: XX bb, XX , 

XX , XX ti and xx W^W~. We can then convert the source significance 

into exclusion curves of the dark matter annihilation cross section for the individual 
dwarf spheroidal galaxies, as detailed in section 4. The data used in this analysis was 
taken while the array was under construction between August 2, 2013 and July 8, 
2014, during which the active detector sites grew from 362 PMTs in 108 WCDs to 
491 PMTs in 134 WCDs. This analysis considered WIMPs annihilating with 100% 
branching ratios into bb, r+r“, ti and W~^W~. 

The 95 % CL limits for the bb and annihilation channels, for the four¬ 

teen dwarf spheroidal galaxies within the HAWC field-of-view, are shown in figure 

Figure shows the exclusion curves for five annihilation channels, for the dwarf 
spheroidal galaxies that HAWC is most sensitive to: Segue 1, Draco, Sextans, Coma 
Berenices and Bootesl. A combined limit was done for each annihilation channel by 
doing a stacked analysis for the fourteen dwarf spheroidal galaxies, and by optimizing 
the total expected counts from all sources. The combined limits are presented both 
in hgures and The HAWC expected fourteen source combined limit was also 
calculated for each annihilation channel and is shown in figure with the hatched 
areas showing the systematic uncertainty on the expected signal from simulation as 
discussed in na. 
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Figure 2: Individual dark matter annihilation cross section limits for the hve best 
dwarf spheroidal galaxies within the HAWC £eld-of-view: Segue 1, Draco, Sextans, 
Coma Berenices and Bootesl. The solid black line in the hrst hve hgures shows the 
combined limit for each annihilation channel as a result of a stacked analysis with 
all fourteen dwarf spheroidal galaxies, while the hatched areas show the expected 
fourteen source combined limits with systematic error on the expected signal. The 
lower right hgure shows the hve annihilation channel combined limits together. 
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6 Summary 

In these proceedings we present individual 95% CL limits on the annihilation cross 
section for fourteen dwarf spheroidal galaxies within the HAWC £eld-of-view. A 
combined limit is also shown from a stacked analysis of all dwarf spheroidal galaxies. 
Limits are presented for a range of dark matter masses and for several dark matter 
annihilation channels, resulting from data collected over a 180 day period during 
construction of the HAWC detector. These are the hrst limits on the dark matter 
annihilation cross section presented with HAWC data. Further analysis of this work 
will include treating the dwarf spheroidal galaxies as extended sources, as well as 
producing exclusion curves for data collected with the entire HAWC array, where the 
limits are expected to improve by about an order of magnitude. 
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